Experimental realisation of quantum illumination 
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We present the first experimental realisation of the quantum illumination protocol proposed in 
Ref.s [S. Lloyd, Science 321, 1463 (2008); S. Tan et al, Phys. Rev. Lett. 101, 253601 (2008)], 
achieved in a simple feasible experimental scheme based on photon-number correlations. A main 
achievement of our result is the demonstration of a strong robustness of the quantum protocol to 
noise and losses, that challenges some widespread wisdom about quantum technologies. 



PACS numbers: 42.50.Ar, 42.50.Dv, 42.50.Lc,03.65.Ud 

Properties of quantum states have disclosed the possi- 
bility of realizing tasks beyond classical limits, originat- 
ing a field collectively christened quantum technologies 
[IHZ]. Among them, quantum metrology and imaging 
aim to improve the sensitivity and/or resolution of mea- 
surements exploiting non-classical features, in particular 
non classical correlations [8HT2 7 ] . However, in most of the 
realistic scenarios, losses and noise are known to nullify 
the advantage of adopting quantum strategies [13 . Here, 
we present the first experimental realization of a quan- 
tum enhanced scheme Q31 [TS] , designed to target detec- 
tion in a noisy environment, preserving a strong advan- 
tage over the classical counterparts even in presence of 
large amount of noise and losses. This work, inspired by 
theoretical ideas elaborated in [T4HT7] (see also[18j). has 
been implemented exploiting only photon-number corre- 
lations in twin beams and, for its accessibility, it can find 
widespread use. Even more important, it paves the way 
to real application of quantum technologies by challeng- 
ing the common believe that they are limited by their 
fragility to noise and losses. 

Our scheme for target detection is inspired by the 
"Quantum Illumination" (QI) idea [21 [T5], where the 
correlation between two beams of a bipartite non- 
classical state of light is used to detect the target hidden 
in a noisy thermal background, which is partially reflect- 
ing one of the beam. In [HI [16] it was shown that for QI 
realized by twin beams, like the ones produced by para- 
metric down conversion, there exists in principle an opti- 
mal reception strategy offering a significant performance 
gain respect to any classical strategy. Unfortunately, this 
quantum optimal receiver, is not yet devised, and even 
the theoretical proposal of sub-optimal quantum receiver 
[19) was very challenging from an experimental point of 
view, and has not been realized yet. 

Our aim is to lead the QI idea to an experimen- 
tal demonstration in a realistic scenario. Therefore, in 
our realization we consider realistic a-priori unknown 
background, and a reception strategy based on photon- 
counting detection and second-order correlation measure- 



ments. We demonstrate that the quantum protocol per- 
forms astonishingly better than its classical counterpart 
based on classically-correlated light at any background 
noise level. More in detail, we compare quantum illu- 
mination, specifically twin beams (TWB), with classi- 
cal illumination (CI) based on correlated thermal beams 
(THB), that turns out to be the best possible classical 
strategy in this detection framework. 

On the one hand our approach, based on a specific and 
affordable detection strategy in the context of the current 
technology, can not aim to achieve the optimal target- 
detection bounds of Ref. [15], based on quantum Cher- 
noff bound [20H!?!?] . On the other hand, it maintains most 
of the appealing features of the original idea, like a huge 
quantum enhancement and a robustness against noise, 
paving the way to future practical application because 
of the accessible measurement technique. Our study also 
provides a significant example of ancilla- assisted quan- 
tum protocol, besides the few previous realizations, e.g. 

In our set up (see Fig. [T]) Parametric Down Conversion 
(PDC) is exploited to generate two correlated light emis- 
sions with average number of PDC photons per spatio- 
temporal mode \i = 0.075, that are then addressed to a 
high quantum efficiency CCD camera (See Supplemen- 
tary Information, Sec. I). In the QI protocol (Fig. [TJi) 
one beam is directly detected, while a target object (a 
50:50 beam splitter) is posed on the path of the other 
one, where it is superimposed with a thermal background 
produced by scattering a laser beam on an Arecchi's ro- 
tating ground glass. When the object is removed, only 
the background reaches the detector. The CCD camera 
detects, on different areas, both the optical paths. In 
the CI protocol (Fig.^), the TWB are substituted with 
classical correlated beams, obtained by splitting a single 
arm of PDC, that is a multi-thermal beam, and by ad- 
justing the pump intensity to ensure the same intensity, 
time and spatial coherence properties for the quantum 
and the classical sources. 

We measured the correlation in the photon number 
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FIG. 1: Experimental setup, a) Quantum illumination: af- 
ter the BBO crystal, where TWB are generated, one of them 
(the "ancilla" ) is reflected towards the detection system. The 
correlated beam is partially detected, together with the ther- 
mal field from the Arecchi's disk, when the object (actually 
a beam splitter) is present, otherwise it is lost (not showed). 
Low-pass filter (95 % of transmission at 710nm) and UV- 
reflecting mirror, not showed, are used to minimize the back- 
ground noise while maintaining low losses. The lens, placed 
at the focal length from the crystal and the CCD camera, re- 
alizes the Fourier transform of the field at the output face of 
the crystal, b) Classical illumination: one beam from PDC 
is stopped and the other one is split at a beam splitter (BS) 
for generating correlated multi-thermal beams. The power 
of the pump is adjusted to obtain the same energy of the 
TWB. c) Detected TWB, in the presence of the object, with- 
out thermal background. The region of interest is selected by 
an interference filter centered around the degeneracy wave- 
length (710 nm) and bandwidth of 10 nm. After selection the 
filter is removed, d) Detected field for split thermal beams in 
the presence of the object, without thermal noise, e) A typi- 
cal frame used for the measurement, where a strong thermal 
background has been added on the object branch. The colour 
scales on the right correspond to the number of photons per 
pixel. 



Ni and AT 2 detected by pairs of pixels intercepting cor- 
related modes of beam "1" and "2" respectively (Fig. [I] 
c-d-e), [26l [28] . With our experimental setup, this corre- 
lation can be evaluated with a single image by averaging 
over all the N pix pixels pairs. Albeit the usage of spatial 
statistic is not strictly necessary, it is practically effective 
and allows to reduce the measurement time (less images 
needed) [27]. 

In order to quantify the quantum resources ex- 
ploited by our QI strategy we introduce a suitable non- 
classicality parameter: the generalized Cauchy-Schwarz 
parameter e = (: SN X SN 3 :> / y/(: S 2 Nx :) (: S 2 N 2 :), 
where (: :) is the normally ordered quantum expectation 
value and S 2 Ni = (Ni — (Ni)) 2 the fluctuation of the pho- 
ton number JVj, i — 1,2. This parameter is interesting 
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FIG. 2: Generalized Cauchy-Schwarz parameter e in the case 
of twin beams, e' TW ', and of the correlated thermal beams, 
e' TH ' , as a function of the average number of background pho- 
tons Nb for a number of background modes M& = 57 (black se- 
ries) and Mb = 1300 (red). The lines represent the theoretical 
prediction at /i = 0.075 (the last estimated independently). 



since it does not depends on the losses and it quantifies 
non-classicality: e < 1 for classical state of light (with 
positive Glauber-Sudarshan P- function), while quantum 
state with negative/singular P- function can violate this 
bound 29: . In Fig. [2] we report the measured e and the 
theoretical prediction. One observes that for TWB e^ 1 ^ 
is actually in the quantum regime (e^ 1 ^ > 1) for small 
values of the thermal background (Nb); in absence of it 
({Nb) — 0) we obtain e^ 1 " 1 ~ 10. As soon as the con- 
tribution of the background to the fluctuation of iVjj be- 
comes dominant, e^ 7 ) decreases quite fast, well below 
the classical threshold. As expected, for THB e^ CI ^ is 
always in the classical regime, and it is equal to one for 
(N b ) = 0. 

We consider an a-priori unknown background, meaning 
that it is impossible to establish a reference threshold of 
photo-counts (usually the mean value of the background) 
to be compared with the possible additional mean photo- 
counts coming from the reflected probe beam (if the tar- 
get is present). Therefore, the estimation of the first or- 
der (mean values) of the photo-counts distribution, typ- 
ical of other protocols (e.g. [H El [H] ) , is here not in- 
formative regarding the presence/absence of the object. 
We underline that this unknown-background hypothe- 
sis accounts for a "realistic" scenario where background 
properties can randomly change and drift with time and 
space. 

For this reason we propose to discriminate the pres- 
ence/absence of the object by distinguishing between the 
two corresponding values of the covariance A12, evalu- 
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FIG. 3: Signal-to-noise ratio (SNR) versus the number of 
background photons (N b ) normalized by the square root of 
number of realization. The red (black) markers refer to 
Mb — 1300 (Mb = 57) and the solid (dashed) theoretical 
curve corresponds to quantum (classical) illuminating beams. 
The estimation of quantum mean values of Eq. |6]) is obtained 
by performing averages of A^ ^ out ^ over a set of JVi mg ac- 
quired images (JVi mg = 2000, 4000, 6000 for twin beams at 
M b = 1300, M b = 57 and thermal beams at M b = 1300, 
respectively). The lowest curve of the classical protocol has 
not been compared with the experimental data because the 
SNR is so low that a very large number of images (out of 
the possibility of the actual setup) is required to get reliable 
points. 



of 



1.000 
0.500 

0.100 
0.050 

0.010 
0.005 

0.001 



(A/ imq =10)- 




20000 



FIG. 4: Error probability P orr of the target detection versus 
the total number of photons of the thermal background Nb 
evaluated with JVimg = 10 (JVi mg = 100 in the inset). The 
black squares and red circles are the data for QI with Mb = 57 
and Mb = 1300, respectively, while red diamonds referes to 
the data for the CI with Mb = 1300. The curves are the 
corresponding theoretical predictions. 



E[X] — ^ Ylk=i represents the average over the set 
of fC realizations that in our experiment correspond to the 
number JC = N p i x correlated pixels pairs. The signal- 
to-noise ratio can be defined as the ratio of the mean 
"contrast" to its standard deviation (mean fluctuation): 



/snr = 
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where "in" and "out" refer to the presence and absence 
of the object. 

For JC » 1, the "contrast" at the numerator of Eq. 
([2]) corresponds to the quantum expected value of the 



covariance i.e. 

. (out) 



<Ai' 



(in)\ 



(SN1SN2), where obviously 



(A\ 2 ) = 0. For a generic prominent background with 
a mean square fluctuation (5 2 Ni,), the "noise" at the de- 
nominator depends only on the local statistical properties 
of the beam 1 and of the uncorrelated background, i.e. 
((5 2 Ai i2 ) ^ (5 2 Ni}(5 2 N b ) (Supplementary Information, 
Sec.II-a) . This is shown in Fig. [5J where the estimated 
covariance of Eq. ([!]) is plotted versus the intensity of 
the thermal background, used in our experiment. As ex- 
pected, the average value of covariance does not depend 
on the quantity of environmental noise, while the uncer- 
tainty bars do. 

While the signal-to-noise ratio unavoidably decreases 
with the added noise for both QI and CI, the quantum 

enhancement parameter (R — /snr //snr) m ^ ne P res ~ 
ence of dominant background and equal local resources 
becomes 



Rr* (6N 1 SN 2 ) QI /(SN 1 SN 2 )ci- 



(3) 



Being R expressed as a ratio of covariances, it is remark- 
ably independent on the amount of losses, noise and re- 
flectivity of the object. 

According to Eq. (JtI the enhancement is lower 
bounded by the amount of violation of the Cauchy- 
Schwarz inequality for the quantum state considered in 
the absence of background, i.e. R ~ s^ 1 ^ /e^ CI ^ > e^ 1 ^ 1 . 
The equality holds for classical states saturating the 
classical bound, Eq = 1 (Supplementary Informa- 
tion, Sec. B). 

In particular, in our experiment we compared the per- 
formance of TWB with a classically correlated state with 
(c I) 

£o =l (hence representing the best possible classi- 
cal strategy), i.e. with a split thermal beams present- 
ing the same local behaviour of the TWB. In this case 
the enhancement can be explicitly calculated obtaining 
R ~ (1 + fJ-)/fJ-, hence the quantum strategy performs 
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FIG. 5: Covariance in the presence, A^j (blue), or ab- 
sence, (green), of the target, a) and b) refer to Ql 
and CI, respectively, for the same number of background 
modes Mb = 1300; c) refers to Ql with a a lower number 
of modes, Mt — 57. Uncertainty bars represent the effect of 
the background noise on the covariance estimation (obtained 
averaging over the iVi mg = 2000, 6000 and 4000 images in 
the graphs from the top to bottom, respectively). Horizontal 
lines are the theoretical values (A^"^ 011 ''}, while the dashed 
lines are the uncertainty interval evaluated theoretically as 



orders of magnitude better than the classical analogous 
when fj, <C 1, namely when a low intensity probe is used. 

Incidentally, since covariance is always zero (i.e. e = 0) 
when using split coherent beams, they do not provide 
a valuable alternative in the situation considered here, 
i.e. when first order momenta are not informative due to 
unknown fluctuating background. 

In Fig. [ij the theoretical prediction for /snr/v^ 
is compared with the experimental data. In perfect 
agreement with theory, the quantum enhancement is al- 
most constant (R > 10) regardless the value of (Nb). 
Therefore, the measurement time, i.e. the number of 
repetitions N mg needed for discriminating the pres- 
ence/absence of the target, is dramatically reduced in 
Ql (for instance, to achieve /snr = 1, ^img is almost 100 
times smaller when quantum correlations are exploited). 

Another figure of merit that highlights the superiority 
of the quantum strategy versus the classical one is the 
the error probability in the discrimination of the pres- 
ence/absence of the target (P err ). In Fig. [4] we report 
P orr versus the number of photons of the thermal back- 
ground (Nb) ■ P C rr is estimated fixing the threshold value 
of the covariance that minimizes the error probability it- 
self. Fig. [4] shows a remarkable agreement between the 
theoretical predictions (lines) and the experimental data 
(symbols), both for Ql and CI strategy. P eII of Ql is 
several orders of magnitude below the CI one and, in 



terms of background photons, the same value of the er- 
ror probability is reached for a value of Nb at least 10 
times smaller than in the Ql case. 

In conclusion, we demonstrated experimentally quan- 
tum enhancement in detecting a target in a thermal ra- 
diation background. Our system shows quantum correla- 
tion (e(QI) ~ 10) even in presence of losses introduced 
by a partially reflective target. Remarkably, even after 
the transition to the classical regime 

( £ (Qi) < 1) due to 

the presence of background ((Nb) 3> 1), the scheme pre- 
serves the same strong advantage with respect to the best 
classical counterpart based on classically correlated ther- 
mal beams. Furthermore, the results are general and do 
not depend on the specific properties of the background 
used in the experiment. 

In paradigmatic quantum enhanced schemes, often 
based on the experimental estimation of the first mo- 
menta of the photon number distribution, such as quan- 
tum imaging protocol [TT], detection of small beam dis- 
placement [5] and phase estimation by interferomctry 
|12j . it is well known that losses and noise rapidly de- 
crease the advantage of using quantum light 13 . This 
enforced inside the generic scientific community the com- 
mon belief that the advantages of entangled and quantum 
state are hardly applicable in a real context, and they will 
remain limited to experiments in highly controlled labo- 
ratories, and/or to mere academic discussions. Our work 
breaks this belief showing orders of magnitude improve- 
ments compared to CI protocol, independently on the 
amount of noise and losses using devices available nowa- 
days. In summary, we believe that the photon counting 
based Ql protocol, for its robustness to noise and losses 
has a huge potentiality to promote the usage of quantum 
correlated light in real environment. 
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ADDITIONAL INFORMATION ON 
EXPERIMENTAL SET-UP 

In our set up a pump laser beam at 355 nm (triplicated 
Nd-Yag Q-switched laser) with 5 ns pulses at 10 Hz repe- 
tition rate, pumps a type II 7 mm long BBO crystal, cut 
for producing collinear degenerate emission. 

The traveling wave PDC generates a spatially mul- 
timode emission, where each mode corresponds to the 
transverse component of a specific wavevector. Pairs of 
correlated modes, corresponding to opposite transverse 
components of the wavevector with respect to the pump 
direction, are found in symmetric positions in the far 
field around the degeneracy wavelength of 710 nm[26- 
28]. Thus, we choose (Fig. 1 c-d-e) two correlated re- 
gions of interests on the CCD array, a Princeton Instru- 
ment Pixis 400 BRW (1340A400 array with pixel size of 
20 ^m, high quantum efficiency, more than 80% and low 
noise in the read out process, about 4 electrons /pixel). 
The proper sizing of the pixels and the centering of the 2- 



dimensional array with sub-mode precision, allows max- 
imizing for each pair of pixels the collection of the corre- 
lated photons and at the same time minimizing the pos- 
sible presence of uncorrelated ones [27,28]. The photon- 
numbcr-correlation, even at the quantum level for QI, 
is realized independently for each pairs of symmetrical 
(translated) pixels belonging to the correlated regions of 
interests of the TWB (THB). The number of pixel pairs 
we used is N pix = 80 with area A p ^ = (480^im) 2 . We 
underline that spatial statistics for the photon-counting- 
based QI protocol is completely unnecessary. Since QI 
does not aim an image reconstruction, two conventional 
photo-counters (i.e. without spatial resolution) will be 
able to accomplish the evaluation of correlation by tem- 
poral statistics, and consequently the QI protocol as well. 
In our proof-of-principle experiment spatial resolving de- 
tectors were employed to have a better control on the 
statistical properties of the field employed, in order to 
provide a proper comparison between the experimental 
data and the theoretical prediction. 

The mean number of photons detected per pixel is 
(N) ~ 4200. The number of collected spatio-temporal 
modes is estimated to be M = 9 • 10 4 by fitting a 
multi-thermal statistics. Thus, taking into account the 
overall transmission and detection efficiency (about 62% 
[27]), the average number of PDC photons per mode is 
/i = (N)/(i]M) ~ 0.075. We measured separately the 
size of the spatial mode, as the FWHM of the correlation 
function between the two beams, A corr = (120 ± 20yum) 2 . 
Therefore, the number of spatial modes is about M sp = 
A p i x /A corr = 16 ± 5 and the number of temporal modes 
Aft = M/M sp = (6 ±2) • 10 3 , the last one being consistent 
with the ratio between the pump pulse duration and the 
expected PDC coherence time, i.e. 1 ps. 

The background field is produced by a pulsed laser 
scattered by a rotating ground glass disk (Arecchi's disk) 
and collimating optics. The spatial properties at the 
CCD plane are set to similar value of the PDC emission, 
while the temporal modes Mf, are selected by adjusting 
the pulse duration and the speed of the disk. 



SNR IN CASE OF PREPONDERANT 
BACKGROUND 

Evaluation of /snr 

For a large number of samples JC >> 1, the "contrast" 
at the numerator of Eq. (2) of the main text corresponds 
to the quantum expected value of the covariance, i.e. 
(A^) - (<Wi<W 2 ), while the mean square fluctuation 
of the covariance (c^Ai^) at the denominator can be cal- 
culated as 
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/C(5 2 A 1)2 ) ~ (S 2 (8N 1 8N 2 )) = (^(5N 1 5N 2 ) 2 ^ — (8N 1 SN 2 ) 2 

( 4 ) 

By replacing SN 2 i-> 5N { 2 " l) + SN b where N^ n) is the 
number of detected photons that has been reflected by 
the target, and N b is the uncorrelated background, the 
right hand side of Eq. Q can be rewritten as 



JC(S 2 A h2 ) 



(5Ni5N 2 ) 2 ) - (8N 1 8N 2 f + (S 2 ^) (5 2 N b ) 
S 2 (SN 1 SNt ) )) + (S 2 ^) (S 2 N b ) , 



where we used the statistical independence of N b and 
the fact that (SNb) =0. It is clear that in the absence of 
the target (situation labeled with the superscript "out"), 
iV 2 (in) = 0, thus (6 2 A < ^ t) ) = (5 2 Ni) (S 2 N b ), since noth- 
ing is reflected to the detector. However, if the the 
background fluctuations (S 2 N b ) is the largest contribu- 
tion to the noise, also when the target is present (indi- 
cated with superscript "in") we can write (8 2 /S^^) ~ 
(5 2 Ni)(S 2 N b y Under this assumption representing a 
realistic situation of a very noisy environment, the SNR 
becomes 



SNR 



(5N 1 5N 2 



(G) 



We underline that Eq. Q holds for a dominant back- 
ground, irrespective of its statistics (e.g. multi-thermal 
or Poissonian). 

In our experiment we consider background with multi- 
thermal statistics. For a generic multi-thermal statistics 
with number of spatiotemporal modes M, mean photon 
number number per mode fi, the total number of detected 
photons is (./V) = Mrjfj, and the mean squared fluctuation 
is (S 2 N) = M7?m(1 + w0 = (N) (1 + (N) /M) [see for ex- 
ample: L. Mandel, E. Wolf, Optical Coherence and Quan- 
tum Optics (Cambridge University Press, 1995)], where 
77 is the detection efficiency 

Thus, the amount of noise introduced by the back- 
ground can be increased by boosting its total number of 
photons (N b ) or by varying the number of modes M b , as 
highlighted from the behaviour of the SNR in Fig. 2, 3 and 
4 of the main text. 



Moreover, both TWB and correlated THB present lo- 
cally the same multi-thermal statistics, but with a num- 
ber of spatiotemporal modes M = 9T0 4 much larger than 
the one used for the background beam (M b = 57 in one 
case and, M b = 1.3 TO 3 in the other). This contributes to 
make the condition of preponderant background effective 
in our realization, even for a relatively small value of N b . 

However, we point out that all the theoretical curves 
reported in all the Figures are evaluated by the exact 
analytical calculation of the four order (in the number of 
^hjo^ns/jgauantum expectation values appearing on the 
rigiit hand side of Eq. @, even if the whole expressions 
are far more complex than the ones obtained with the 
assumption of preponderant background. 



Quantum enhancement R 

Starting from Eq. ([6| and considering the same 
local resources for classical and quantum illumination 
beams (in particular the same local variance (<5 2 7Vi) CI = 



Qi 



1,2)) the enhancement of the quantum 



protocol can be easily obtained as 
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with e = (: 5NxSN 2 :) / y/(: <5 2 iVi :)(: S 2 N 2 :) being 
the generalized Cauchy-Schwarz parameter introduced 
in the main text. The covariance of two correlated 
beams obtained by splitting a single thermal beam is 
(6N 1 6N 2 )th = Mj71r72Atl.fr, whi le the one of TWB is 
(SNi8N 2 )tw — Mr)irj2[TTw(l + firw) (see for example 
Rcf. [27]). By using this relation with the assumption 
of the same local resources, hth = Htw = (i we can 
derive explicitly R ~ (1 + which is insensitive to 

the amount of noise and loss. On the other side the gen- 
eralized Cauchy-Schwarz parameter for a split thermal 



beam is e 



(CI) 



1, where the subscript "0" stands for "in 



absence of background" , as it can be easily derived from 
the equations of covariance and single beam fluctuations 
used previously. 



As described in the text e 



(CI) 



1 represents the best 



result for classical states and, in this sense, the compari- 
son with split thermal beams represents the comparison 
with the "best" classical case. 



